Necessary Condition for the Quantum Adiabatic Approximation 
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A gapped quantum system that is adiabatically perturbed remains approximately in its eigenstate 
after the evolution. We prove that, for constant gap, general quantum processes that approximately 
prepare the final eigenstate require a minimum time proportional to the ratio of the length of the 
eigenstate path to the gap. Thus, no rigorous adiabatic condition can yield a smaller cost. We also 
give a necessary condition for the adiabatic approximation that depends on local properties of the 
path, which is appropriate when the gap varies. 



The quantum adiabatic theorem asserts that a contin- 
uously perturbed and gapped quantum system remains 
in its instantaneous eigenstate in the limit where the rate 
of change of the perturbation vanishes [T] . This assertion 
is quantified via the adiabatic approximation [2], which 
provides a relation between the rate of change of the per- 
turbation and the fidelity of the evolved state with the 
final eigenstate. The adiabatic: approximation is a key 
part of quantum computing as it determines the com- 
plexity of several quantum algorithms [3j |4] |5j [6]. In 
fact, adiabatic quantum computation [3 ], in which the 
result of a problem is encoded in the ground state of a 
(final) Hamiltonian, is equivalent to standard quantum 
computation [7]. Further, adiabatic approximations play 
an important role in areas like Born-Oppenheimer theory, 
the quantum Hall effect and STIRAP [5]. 

It is important to remark that some familiar adia- 
batic approximations are known to be insufficient |Hj and 
sometimes unnecessary [101 111] . The growing interest on 
the adiabatic approximation has spurred work on corre- 
sponding rigorous conditions jTUJ [TJ] . In this manuscript 
we give a rigorous lower bound for the evolution time (or 
cost) of adiabatic processes that prepare the final state. 
This bound is also valid for more general quantum evo- 
lutions [5J E DUES. 

Let {H(r)}, with r £ [0,1], be a given continuous 
Hamiltonian path and {IVK 7 *))} the corresponding non- 
degenerate eigenstate path (eigenpath). Adiabatic evo- 
lutions aim to prepare IV'(l)) a t bounded precision from 
IV'(O)) by choosing a proper schedule r(t). We recently 
argued |llj that the relevant quantities for the adiabatic 
approximation are not only the minimum eigenvalue gap 
of the Hamiltonians, A, but also the length of the path 
to be traversed, L. We presented a method that adi- 
abatically prepares the final state by evolving with the 
Hamiltonians for suitable random times [TT]. The aver- 
age cost of the randomization method is 0(L 2 / A) when 
the rate of change of the eigenstates along the path and 
the corresponding eigenvalues are known. A more effi- 



cient method to traverse the eigenpath for this case was 
introduced in [B]. This method uses Grover's fixed point 
search and the schedule r(t) is non-monotonic. It results 
in a cost (D(L(log L) 2 /A). (We use the soft order nota- 
tion O to hide doubly logarithmic factors.) Finally, we 
recently derived a non-monotonic quantum algorithm or 
process that dynamically estimates the rate of change of 
the eigenstates, and results in a cost C(LlogL/A) under 
broader assumptions |13j . 

Here we also consider general quantum processes that 
prepare the final eigenstate from the initial one, at 
bounded precision, by evolving with the Hamiltoni- 
ans. Thus, we do not exploit the unknown structure of 
{H(r)}, rather we work in the so-called black-box model 
where the only assumption is to be able to evolve with 
H(r(t)) for some schedule r(i). We then prove that the 
cost of such processes is, at least, 0(L/A). 

To prove such a lower bound on the cost we introduce 
particular instances of Hamiltonian paths {H(r)}, and 
reduce them to problems for which a query-complexity 
bound is known or can be easily obtained. For example, 
to show the scaling with A, we can simply consider the 
adiabatic version of Grover's search EJ El [TS]. If 
N is the problem size, the minimum gap is \/\N and 
L < tt in this case. The lower bound y/N G C(l/A) 
for the search problem is a celebrated result in quantum 
computation |16j . Showing the dependence of the cost on 
the path length and the minimum gap requires a different 
analysis that constitutes our main contribution. 

For the instances considered below, the relevant 
Hamiltonian-cigcnvalucs, gaps, and rates of change of the 
eigenstates are known and remain constant for all r. We 
also clarify that a better bound may be obtained if these 
quantities vary along the path. Still, we show that any 
rigorous adiabatic approximation based on local prop- 
erties of the path cannot yield a schedule that satisfies 
r(t) > cA(r)/\\ \d r ip(r))\\ for all r, and specific c > given 
below. 
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Before obtaining the necessary condition for the adia- 
batic approximation we give a precise definition of L and 
comment on the resulting cost for the worst-case. In the 
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instances considered below the path length is 

L= A \dMr)) \\dr. (1) 
Jo 

(See [TT] for a general definition of L.) With no loss 
of generality we assume (d r tp(r)\ip(r)) =0. L is the 
only natural length in projective Hilbert space (up to ir- 
relevant normalization factors). An upper bound on L 
is H-H1I/A, with \\H\\ = max ||<9 rJ ff||, and ||.|| the oper- 
ator norm. In the worst-case this bound is tight and 
L/A ~ ||9 r if ||/A 2 . However, in many cases of interest L 
can be bounded independently of A and the algorithms 
in [6l HH [13] result in much smaller implementation costs 
than those determined by other rigorous adiabatic ap- 
proximations [TUJ Q2] . 

The remainder of this manuscript is organized as fol- 
lows. We first give a simple proof of the necessary con- 
dition for the adiabatic approximation when the eigen- 
states are degenerate. This proof will motivate the de- 
sired result in the non-degenerate case, proven later. 

Eigenpath within a degenerate subspace. We 
use constant-gap piece- wise adiabatic universal computa- 
tion |17j . Let U = U n ■ ■ ■ Ui be a unitary quantum circuit 
acting on a Hilbert space H, initially in For I > 1, 
we define the parametrized Hamiltonians (see Fig. [IJ 

Hi{s) = ~ [C0fl(7r*)(l ® (\l - 1) {I - 1| - |0 (Z|) (2) 

+ sm(irs)(Ui <g> |/) (I - 1\ + U\ <g> \l - 1) (Z|)] , 

acting on 7i ® X, with s 6 [0,1] and A the relevant 
gap. The Hilbert space X encodes the step of the cir- 
cuit using the so-called clock states [7j [18] . For arbitrary 
we define | 7i ) = (U l ---U 1 10 and | 7o ) = |^) |0). 
Each i?;(s) leaves invariant the subspace spanned by 
{N-i) > | 7 z)}- The eigenstates of Hi(s) are degenerate; 
however, we are interested in the continuous path de- 
termined by the eigenstates \<f> (s)) — cos(7rs/2) | 7 i-i) + 
sin(7rs/2) |7j). 

The continuous Hamiltonian path {H(r)}o< r <i is 
constructed by concatenation of n path segments, 
each segment linking | 7 _i) to | 7 /). Then | -0 ( 1 ) ) = 
(U n ■ ■ ■ U± \n) and the final eigenstate of the path 
contains the state prepared by the circuit. The I— th 
path segment, where (I — l)/n < r < l/n, is determined 
by H(r) = Hi (rn — I + 1) (Fig. tt). The relevant eigen- 
state for this path is \ip(r)) — \<p r {rn — I + 1)). Because 
the rate of change of the eigenstate is constant in the 
intervals (I — X)/n < r < l/n, and each path segment 
occurs in a two-dimensional subspace, the path length is 
L = im. The adiabatic procedure considered is universal 
for quantum computation. 

No general unitary quantum process that evolves with 
the Hamiltonians {H(r)} for some schedule r(t), and in- 
terleaves these evolutions with other (known) operations, 
can prepare the final state at bounded precision with cost 
less than 0(L/A) g 0(n/A). This lower bound easily 
follows from considering those quantum circuits U built 



Hi(p) #i(i) = #m(o) 
c 

|7i-i) hi) = Pibl-i) 

FIG. 1: A quantum process that traverses the I— th path seg- 
ment applies Ui. 

from particular instances of Grover's search. Each op- 
eration Uj that composes the circuit is then a so-called 
Grover iteration |15j which uses one search query. The 
size of the circuit is n — y/N, where N is the size of 
the problem. In [TS] it was shown that the cost of any 
continuous-time quantum algorithm that uses the Hamil- 
tonians H(r) and outputs the desired state is, at least, 
T G 0(VN I A); otherwise the search problem could be 
solved using less than \[N queries. Since 
this proves the result. A basis for the eigenspace of the fi- 
nal Hamiltonian is {\x} \n)} x , and thus reaching the final 
eigenspace is easier than drawing up the right state. 

Non-degenerate eigenpath. We now consider the 
case where the eigenstates are non-degenerate. We show 
that, for any given {H(r)} and A, no generic evolution 
induced by a Hamiltonian H(r(t)) + Ho{t), where ffo(i) 
is a driving Hamiltonian, can prepare an approximation 
to the final eigenstate in time less than 0(L/A). The 
result is valid for any r(t) and H]j(t) that do not depend 
on the unknown structure of H(r), but may depend on 
the known gap and rate of change of the eigenstates [20] . 
This information is usually not available when deriving 
adiabatic approximations; our lower bound on the cost 
clearly encompasses those cases as well. Our setup is gen- 
eral, as it may include controllcd-Hamiltonian evolutions 
as well as intermediate measurements. To prove the re- 
sult we consider instances that reduce to instances of the 
ordered search problem [3T] . A discrete-query complexity 
lower bound, based on the adversary method [16, 22, 23 , 
is known for this problem. We extend this result to the 
continuous-time query setting. 

Let x — [xq, . ■ ■ , x n -i], Xj <G {0,1}, be a secret word 
(input) and define the Hamiltonians 

H l x = A\x(l), +,..., +)(x(l), +,..., + | , (3) 

acting on n qubits, where x(l) — [xq, . . . , xi-i], and 
|+) = (|0) + |1))/V2. The evolution induced by H l x for 
time it/ A is equivalent to a phase query Q l x that, for 
input state |ao, • • • , a/-i) , outputs — |ao, ■ • • , ai-i) when 
cij = Xj V < j < I — 1, or does nothing otherwise. Basi- 
cally, with one query Q x we answer the question: are the 
first I bits of x equal to [clq, . . . , a/_i]?. The Hamiltonian 
path we will consider is one that has \x(l), +,...,+) as 
intermediate eigenstates. Then L £ 0(n). As in the de- 
generate case, this path is built upon n segments, each 
interpolating H l x with H l x +1 (see below). A measurement 
on the final state allows us to learn the input. 

We now reduce our problem to the ordered search prob- 
lem. The goal of ordered search is to find the input 
x that represents a marked item in an ordered list of 
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N = 2 n elements. A query in this case, acting on input 
a = [ao, . . . , a n _i], indicates whether the marked element 
x is before position a or not. The corresponding phase 
query R x puts a ± sign on the input state depending on 
the answer (see below). It requires 0(n) queries of type 
R x to solve the problem [22] . We remark that a query of 
type Q l x can be implemented using two queries of type 
R x : to decide if the first / bits of the secret word coincide 
with a(l) = [ao, . . . , oj-i], it suffices to call R x with in- 
puts [ao, . . . , aj-i, 0, . . . , 0] — 1 (binary subtraction) and 
[a , . . . , di-i, 1, . . . , 1], respectively. If the corresponding 
R x 's outputs are different, then a(/) = x(l). Then, 

Q l x = U l R x V l R x W l , (4) 

where U , V , and W l are the ^-independent unitaries 
used to build a circuit that simulates Q l x using R x . 

A simple argument based on the above analysis intu- 
itively explains our result. It is now clear that no gen- 
eral quantum process that uses type-Q^ queries can find 
the secret word with less than 0(n) queries; otherwise, 
by replacing the queries, we would solve ordered search 
with less than 0(n) type-R x queries. Since each type- 
Q l x query can be implemented using the Hamiltonians 
H l x for time 7r/A, it is plausible that the total cost of a 
continuous-time quantum process that uses H x is at least 
0(n/A) e 0(L/A). In fact, using the equivalence be- 
tween continuous- and discrete-time query models in |24j 
yields a lower bound for the cost 0[L/(A log(L/A))]. 
Our formal proof below, that uses a version of the ad- 
versary method in the continuous-time setting, will avoid 
the logarithmic correction in the cost. 

Let H l x (s) = cos(7rs/2)iJ^ 1 +sin(7rs/2)iJ^ . We build 
a particular Hamiltonian path as in the degenerate case: 
for (I - l)/n < r < l/n we set H x (r) = H l x (rn - I + 1). 
The cigenstates \ip(r)) of H x (r), with lowest eigenvalue, 
are non-degenerate. Further, the gap of H x (r) is A for all 
r. A quantum algorithm that traverses the eigenpath will 
aim to prepare \x(l), +, ...,+) from \x(l — 1), +,..., +), 
learning a bit of information about the secret word at 
each path segment. The path length is L = nn/2. 

We consider evolutions with Hamiltonian H x (r(t)) + 
Hn(t), where r(t) and H£>(t) do not depend on x, the 
only unknown quantity in this case. Ideally, after some 
time T > 0, the initial state approximately evolves to the 
desired eigenstate \ip(l)) — \%o, ■ ■ ■ , i„_i). Evolutions 
of this type include continuous-time processes based on 
eigenpath traversal, such as adiabatic evolutions where 
Hq = 0. If \<fi x (t)} is the evolved state, 

id t \4> x {t)) = (H x (r(t)) + H D {t)) \<f> x (t)) . (5) 

The initial state ^(O)) is also independent of x. 

We use the adversary method to show that T is at 
least 0(n/A). At its core, the adversary method pro- 
vides a bound for the rate of change of the overlap be- 
tween evolutions corresponding to different inputs. To 
distinguish between these inputs, the evolved states must 
satisfy \{<p x (T)\(f> y (T))\ < e, for some T > and small e. 



Let r be an adversary matrix: an irreducible symmet- 
ric matrix with non-negative entries and zeros in the di- 
agonal. Denote by ||r|| its operator norm, and by v the 
principal unit eigenvector, Tv = |L||u. The following 
function serves as a measure of the distinguishability be- 
tween evolutions with different inputs: 

(t)\4>y{t)) . (6) 

Since the initial state is independent of the input, W(0) = 
\\T\\. Moreover, W(T) < e\\T\\. 

An upper bound on |9tW(£)| can be obtained if we 
extend known results for the ordered search problem to 
our case. We write the Hamiltonians H l x as (see Eq. Q) 

H l x = A(l - Q l x ) = A(l - U l R x V l R x W l ) . (7) 

Then, using Eq. ^ and the triangle inequality, 

\d t (Mt)\Mt))\ <^\(Mt)\H l x -H l v \^(t))\ (8) 
< 2Amax\{4> x (t)\ U l R x V l R x W l - 
-U l R y V l R y W l \^ y (t)) \ . 

For each x, we define a 2™-bit string a x with en- 
tries a x = 0, if i < x, and a x = 1 otherwise. Let 
Pi = \i){i\ and write {(j) l x {t)\ = (cj) x (t)\U l R x V l , and 
|^(t)) = V l RyW l \<t> y {t)). The triangle inequality and 
Eq. (|| yield 

IW*(*)l0ir(*)>l <2Amax £ [K0L(*)|W \<t> y {t)}\ 

+ \(Ut)\U l Pi\fty(t))\)- 

For i £ {0, . . . , 2™ — 1}, we introduce the matrices T l as 

pi _ j ^x,y if ol x cty , 
x ' y [0 otherwise . 

Then, 

\d t W(t)\ < 2A ri >v «x«vniax[|<^(t)| P t W l |0,(t))| 

i,x,y 

+ \(ut)\u l m l y (t))\} 

< 4Amax||r l || . (10) 

i 

In addition, (1 - e)||r|| < W(0) - W(T) and 

W(0) - W(T) < [ \d t W(t)\dt (11) 
Jo 

< T4Amax||r l || . 

i 

The spectral lower bound for the cost of the pro- 
cess that approximates the final eigenstate is T G 
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O (||r||/(4Amaxj ||r z ||)). We use the adversary matrix 
for the ordered search problem [22] : 



r _ J Hd(a <t 



,a v ) 



if 7^ a y 
otherwise 



(12) 



where Hd(a x ,a y ) is the Hamming distance. This choice 
yields ||r|| > n and maxjr 1 !! < n. Thus, T G 
f2 (n/(47rA)). Since L = irn/2, this proves the result 
TeO(L/A). 

With a similar construction we can prove a lower 
bound on the cost of continuous-time query algorithms 
that solve the ordered search problem, obtaining T £ 
f2(n/(27r)) in this case. 

In the derivation of our result we used the fact that 
the gaps of {H x (r)} are constant along the path. Nev- 
ertheless, the Hamiltonians in the path may have differ- 
ent (known) gaps and one could be interested in design- 
ing algorithms for cigenpath traversal with a schedule 
r(t) that depends on the local gap A(r). Consider again 
the instances above and redefine H x {r) <— q(r)H x (r), 
for known q(r) > 0. The new Hamiltonians have gaps 
A(r) = q(r)A. We can replace A by A(r(t)) in Eq. Q. 
Following the steps above, and using the adversary ma- 
trix of Eq. (Il2|) , we obtain 



(13) 



(1 - e)n < 4tt [ dt A(r(i)) 
Jo 



We prove by contradiction that adiabatic approxi- 
mations based on local properties of the path cannot 
yield a schedule satisfying the local condition r(t) > 
cA(r) /\\\d r ip(r)}\\, for some c > and all r. If such con- 
dition is satisfied, the inverse function t(r) exists and 
Eq. (pi) yields 



(l-e)n < 



47T 



4ttL 



dr|||0 r ^(r)) || = 
o c 



(14) 



Since L = wn/2 in this case, it is clear that c < 2tt 2 / (1— e) 
or otherwise the above inequality is inconsistent. 

Conclusions. We proved that no general quantum 
process that approximately prepares the eigenstate of a 
Hamiltonian H{\), by evolving with the path of Hamil- 
tonians {H(r)}o< r <i, for any schedule r(t), can achieve 
its goal in time less than 0(L/A). The same bound on 
the time applies also for more general evolutions with ad- 
ditional driving Hamiltonians. Interestingly, some quan- 
tum processes for eigenpath traversal almost achieve the 
bound under some assumptions |13j . We also gave a nec- 
essary local condition, valid even when the gaps of the 
Hamiltonians and rates of change of the eigenstates are 
known along the path. In this case we proved that if the 
schedule satisfies r(t) > cA(r) /\\\d r ip(r)}\\, for a specific 
c > and all r, the quantum process will not succeed in 
the state preparation. If only A < min r A(r) is known, 
our result suggests that no general monotonic schedule 
yields a cost better than 0(L/A). 
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